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The event-by-event dynamical fluctuations in kaon-to-proton and proton-to-pion ratios have 
been studied in dependence on center-of-mass energies of nucleon-nucleon collisions -^s. Based 
on changing phase space volume which apparently is the consequence of phase transition from 
hadrons to quark-gluon plasma at large ^/s, the single-particle distribution function / is assumed 
to be rather modified. Varying / and phase space volume are implemented in the grand-canonical 
partition function, especially at y/s > 17 GeV, so that hadron resonance gas model, when taking 
into account the experimental acceptance A and quark phase space occupation factor 7, turns to be 
able to reproduce the dynamical fluctuations in {K^ + K^)/{p + p) and {p + p)/{-K'^ + n~) ratios 
over the entire range of ^/s. 
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I. INTRODUCTION 



In order to characterize the fluctuations and collective properties in heavy-ion collisions, phase structure and 
event-by-event fluctuations have been suggested [ll-|Jl ■ From theoretical point -of-view, the dynamical fluctuations are 
O ■ particularly essential regarding to examining of existing statistical models [^, , characterizing of particle equilibration 
' I [7| and the search for unambiguous signals of new state of matter [8|, |9| . From the phenomenological point -of-view, 
p^, they have been measured at SPS and RHIC energies |10l - [l3 |. i.e, 6 < a/s < 200 GeV. 

(-H ' The success of hadron resonance gas model (HRG) in describing the dynamical fluctuations in kaon-to-pion ratios 

[l| and the ingredient to have other implications are the motivations of this work, where the fluctuations in the 
remaining particle ratios, that have been measured so far, i.e, kaon-to-proton and proton-to-pion, shall be analyzed. 
To bring theoretical and experimental data very close to each others, the experimental acceptance A, grand-canonical 
statistics and quark phase space occupation factor, 7. have been utilized. It has been noticed that the fluctuations 
QQ ■ over the whole range of center-of-mass energies ^/s exhibit a non-monotonic behavior and a certain conflguration of 
l/~j . 7, A and grand-canonical statistics, are not able to reproduce the whole experimental data. 

■TJ" ' The concept of phase space dominance suggested by Fermi [TjI six decades ago is applied in this model. The 

r — , assumption of equilibrium single-particle distribution function / and extensive thermodynamics that have been shown 
f^ ' to perfectly reproduce all thermodynamic quantities, including multiplicities and fluctuations, as long as the energy 
^^ , density is not high enough to derive the hadronic system into the new state of matter, quark-gluon plasma (QGP) 
^-H . should be no longer valid, particularly, when the energy density exceeds the critical value, as the case at RHIC and 
^ I LHC energies. Across the phase transition, the phase space volume, symmetries and accordingly the effective degrees 
of freedom are likely subject of a prompt change. Therefore, the phase space volume, in which the microscopic states 
are distributed according to / function is conjectured to be a subject of modiflcation, as well. In present work, 
changing phase space volume and / are implemented in the grand-canonical partition function, explicitly at y/s > 17 
GeV, from which the particle number and fluctuation have been calculated. 

There are many reasons to apply HRG model. This model provides a good description for the thermodynamical 
evolution of the hadronic system below the critical tenrperature [a,[ij,[ia| and has been successfully used to characterize 
the conditions deriving the chemical freeze-out [1^, Il3 • Apparently, there are many reasons speak for studying of 
dynamical fluctuations in the baryon-meson yield ratios. Statistically, the interplay between fermions and bosons 
is crucial. The energy threshold required to create baryons and mesons is of a great interest. Size of the system 
is another ingredient of this study. Phenomenologically, the strange fluctuations, especially, when passing through 
deconfinement phase transition, are conjectured to affect the dynamical fluctuations in the hadronic final state. The 
strange dynamical fluctuations are expected to survive through the mixed phase. Therefore, these particle ratios are 
sensitive to the symmetry change and the dynamics of deconfinement and chiral phase transitions, respectively. 
The paper is organized as follows. The model is given in section|lTl where non-extensive single-particle equilibrium 
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distributions arc introduced. Section Hill is devoted to the dynamical fluctuations in kaon-to-proton and proton-to- 
pion ratios. Discussion and final conclusions are elaborated in section ITVl 

II. MODEL 

A. Single— Particle Equilibrium Distribution of Hadrons 

The grand-canonical partition function is given by Hamiltonian and baryon number operators, H and h, respectively, 
and depends on temperature T = 1//3 and chemical potential /x, 

Z{P,V,^i) = Tr [7 exp'^(^^~^)j . (1) 

It can be characterized by various but a complete set of microscopic states and therefore the physical properties of 
the quantum systems turn to be accessible i n. app roximation of non-correlated free hadron resonances. Therefore, 
the resonances are treated as a free gas p, [ij, [l^, [l3, [3- Each of them is conjectured to add to the overall 
thermodynamic pressure of the strongly interacting hadronic matter. This statement is valid for free, as well as, 
strong interactions. It has been shown that the thermodynamics of strongly interacting system can be approximated 
to an ideal gas composed of hadron resonances with masses < 2 GeV [18, _20|. The heavier constituents, the smaller 
thermodynamical quantities. The main motivation of using the Hamiltonian is that it contains all relevant degrees of 
freedom of confined and strongly interacting matter. It implicitly includes the interactions that result in the formation 
of new resonances. In addition, this model has been shown to provide a quite satisfactory description of the particle 
production and collective properties in heavy-ion collisions. 

The conservation of baryon number rif, represents an additional constrain on grand-canonical partition function. 
Therefore, Eq. ([T]), can be re-written as 

Zg,{T,V,y) = TV [exp^'^-"! , (2) 



exp 1 y 

where a, besides /3, are Lagrange multipliers in the entropy maximization. The physical meaning of a is controller 
over number of particles in the phase space, i.e, acting as chemical potential. 

a = ln£-lnT-ln7V. (4) 

It combines intensive variables, T and N with an extensive one £ = X]"5i^*j where gi and e^ are degeneracy factor 
and energy of z-th cell in the phase space, respectively. The most probable state density is to be found by Lagrange 
multipliers, where one of them, a, has been expressed in term of the second one, /3, and the occupation numbers of 
the system. Apparently, a gives how the energy £ is distributed in the microstates of the equilibrium system and 
therefore, can be understood as another factor controlling the number of occupied states, at the microcanonical level. 
Under these assumptions, the dynamics of partition function can be calculated as a summation over single-particle 
partition functions Z'^g^ of all hadrons resonances. 

/>oo 
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where ± stands for bosons and fermions, respectively. A^ = exp{pi/T) is the i-th particle fugacity. As given above, 
7 = exp(— a) is the quark phase space occupation factor. 

As obtained in Ref. [21[, the equilibrium distribution function / is conjectured to be no longer valid, especially, 
when the energy density is high enough to derive confined hadronic into deconfined partonic matter. The equilibrium 
is settled, when d-tS = Q and the probability current entirely vanishes, as well (22|, |23|. 

f{x,k,t) ~ f,g{x,k) Q{x,k), (6) 

where Q{x, k) is conjectured to reflect the change in phase-space when the hadronic degrees of freedom are replaced 
by partonic ones. It can be interpreted as a measure for the non-extensivity. Therefore, the partition functions reads. 
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where single-particle distribution function and the particle number can be written as 
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III. DYNAMICAL FLUCTUATIONS IN BARYON-MESON RATIOS 

The fluctuations in particle number, Eq. ([9]), are mainly given by the susceptibility, which is the derivative of 
particle number (n) wrt chemical potential fi. 
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where (• ■ ■ ) stands for statistical average. 

At final state, i.e, when the chemically relaxing system absolves the chemical freeze-out process, the resonances 
are conjectured to decay either to stable particles or to other resonances. This chemical process has to be take into 
account in the particle numbers and fluctuations given above as follows. 
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where hj^i is the branching ratio for the decay of j'-th resonance to i-th particle. In order to characterize the stage 
at which the chemical freeze-out takes place, we use the model introduced in Ref. [la, [13, ^M,- The ratio s/T^, 
where s = S/V is the entropy density, is assigned to a constant value. This thermodynamical condition perfectly 
reproduces the freeze-out line over a wide range of chemical potential, /x, which in turn can be empirically expressed 
in center-of-mass energy, ^/s in GeV, |24| 
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The event-by-event fluctuations in a ratio of two particles ni/n2 are [j| 
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which includes both dynamical and statistical fluctuations. The third term of Eg . [TS] counts for fluctuations from the 
hadron resonances which decay into particle 1 and particle 2, simultaneously. In such a mixing channel, all correlations 
including quantum statistics ones are taken into account. Obviously, this decay channel results in strong correlated 
particles. To extract statistical fluctuation, we apply Poisson scaling in mixed decay channels. Experimentally, there 
are various methods to construct statistical fluctuations J12l | . Frequently used method is the one that measures particle 
ratios from mixing events. 
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Subtracting Eq. [inifroni Eg. 1151 then the dynamical fluctuations in ni/n, ratio read 
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IV. DISCUSSION AND CONCLUSIONS 



As introduced previously, HRG assumes that the hadron resonances are point-hke and non-correlated free gas. 
Therefore, it is conjectured that the average multiplicity (n), Eq. (J12p . and the dynamical fluctuations cr^ ,j^ , 

Eq. (|15p. are not strongly dependent on the volume fluctuations. Therefore, we assume that volume fluctuations over 
the entire range of ^/s are minimum and thus neglected. The experimental acceptances of the different detectors A 
have been taken into account. Also, the quark phase space occupation factor, 7, has been estimated. The dynamical 
fluctuations in particle yield ratios with and without Q are plotted in Figs. [T]&[21 
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FIG. 1: Dynamical fluctuations in (p + |5)/(7r+ + tt^) ratio as function of centcr-of~mass energy yG. Taking into 
account the experimental acceptance A, the lower curve represents HRG results at finite value for 7. It reproduces 
perfectly the SPS data. The dashed region shows that RHIC data are largely underestimated. The upper curve 
represents HRG that reproduces HSD simulations. RQMD and HSD have been proceeded using RHIC configurations. 
When phase space factor Q is switched on, HRG matches with RHIC data. The overall curve is given in the right 
panel. 
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FIG. 2: Dynamical fluctuations in {K^ + K~)/{p + p) ratio as function of center-of~mass energy ^/s. SPS fluctuations 
rapidly decrease to negative values, while RHIC fluctuations are positive and slowly increase with increasing ^/s. 
Taking into account the experimental acceptance A, the lower curve represents HRG results at finite value for 7, that 
satisfactorily match with SPS data. The dashed region shows that RHIC data are largely underestimated. The upper 
curve represents HRG results, when phase space factor Q is switched on. RHIC data perfectly matches with this 
curve. The curves arc not stemming from statistical fitting process. 



Left panels in both figures show a full comparison between different data using different parameters. Right panels 
give the final results that satisfactorily describe the non-monotonic behavior. The quark phase space occupation 
factor 7 alone is apparently not able to reproduce the whole data set, although the experimental acceptances of the 



different detectors have been taken into account. The largest disagreement occurs at ^/s > 17GeV, is given by bottom 
curves, especially the dashed parts. RQMD and HSD are used to generate dynamical fluctuations in corresponding 
particle yield ratios using RHIC configurations. Therefore, their agreement with RHIC data is satisfactory. The 
dynamical fiuctuations in same particle yield ratios and strangeness production in HSD transport approach from SIS 
and SPS energies has been given in Ref. [2^. Although, both RQMD and HSD refiect a global behavior that the 
fluctuations alternate between negative and positive values, RQMD seems to work nicely at very high energies while 
HSD obviously covers a much wide range of energy. 

Some details about RQMD and HSD are now in order. The relativistic quantum molecular dynamics (RQMD) 
is based on color rope mechanism and excitation and/or fragmentation of color strings followed by hadronic re- 
interactions. The color rope mechanism represents an important source for the production of strange hadrons. The 
comparison with experimental data shows that RQMD reproduces very well the main properties in hadron production 
and collective properties. Therefore, it is believed to be a reliable model for generating the full momentum space of 
hadron resonances. The particle yields from RQMD would be accompanied with statistical uncertainties that basically 
can be kept arbitrarily small, i.e, generation of large ensemble. That the microscopic state of RQMD at chemical 
freeze-out is in good agreement with chemical and kinetic equilibrium leads to believe in RQMD chemical freeze-out 
parameters [2a |. 

The HSD transport approach is based on various degrees of freedom, including quark, diquark, string and hadron. 
It describes the creation of dense/hot hadronic matter, the in-medium modification of hadron properties and the 
overall dynamics. HSD gives numerical solution of a coupled set of relativistic transport equations for particles 
with in-medium self-energies of test-particle. It applies two approaches. One at high energy, where the inelastic 
interactions are described by MC-techniques, FRITIOF string model. At low energy, modeling of hadron collisions is 
based on experimental inputs. The transport approach implemented in it takes into account formation and multiple 
re-scattering of leading pre-hadrons and hadrons and has been d esig ned to reproduce nucleon-nucleon, meson-nucleon 
and meson- meson cross section data in a wide kinematic range |25| . 

In addition to the assumptions introduced in [23, [23, we give here a novel one. We assume that the prompt 
raise at ^/s ~ 17 GeV is to be understood according to a modification in the phase space volume. To this end, it 
has been concluded in Ref. [ij that / is a subject of modification, especially at large \/s. In this limit, the energy 
density available to the system turns to be high enough to cause the hadronic matter, where equilibrium / is perfectly 
able to reproduce almost all essential transport properties and thermodynamic quantities, including the dynamical 
fluctuations, to go through a phase transition into QGP. Such a phase transition apparently results in various types 
of modiflcations, such as symmetries and degrees of freedom. Also the conflgurations of microstates in phase space 
volume <Px<Pk and the single-particle distribution function / are not exceptions. 

Here, we apply this model to the dynamical fluctuations in {K'^ + K~)/{p + p) and [p + p) / [t^^ +7r^) ratios. These 
particle yield ratios combine both strangeness and light boson fluctuations. Additionally, they are very sensitive to the 
deconfinement and chiral phase transitions, respectively. RQMD and HSD simulations result is positive fluctuations at 
high energy. Implementing 7 and Q in the grand canonical partition function of HRG results in a very well description 
of the experimentally measured fluctuations in these particle yields over the entire range of -^5. Although the 
experimental fluctuations are negative overall, the non-monotonic behavior is very well reproduced over a wide range 
of \/s. It is obvious, that these dynamical fluctuations non-avoidablely refer to non-extensive and non-equilibrium 
state of matter, that basically differs from the one at SPS energies. The modification of the state of matter has been 
combined in the factor Q{x,k) ~ exp(2a;a;A:) whose numerical value is given in Tab. U and kept unchanged with 
changing y/s, right panel of Figs. [T]and[2] 

Left panel of Fig. [T] includes SPS and RHIC experimental data (circles) with RQMD and HSD simulations (rect- 
angles) for the dynamical fiuctuations in {K^ + K~)/{p + p) ratio. The results from HRG model are given by the 
curves. It is clear that SPS data can be reproduced by HRG using a certain set of 7 and A parameters. The whole set 
of parameters is given in Tab. H] The parameters suitable for SPS are no longer able to reproduce RHIC data. To this 
end, Q has to be switched on. It is clear that HSD simulations can be reproduced in HRG, while RQMD not overall, 
especially in middle region. Almost the same behavior is present in the left panel of Fig. [2] for (p + p)/{tt~^ + 7r~) 
ratio. Here, the fiuctuations are alternating between positive and negative values. Even at low energy, <7dyn decreases 
and flips its positive sign, with increasing ^/s. At y/s > 17 GcV, crdyn suddenly jumps into the positive region. This 
prompt sign-exchange can be taken as an order parameter. When the energy available to the system is high enough 
to create protons, the dynamical fluctuations get positive. The parameters 7 and A that have been used to reproduce 
SPS data turn to underestimate the RHIC results, the dashed curve. Implementing Q produces the top curve, which 
matches satisfactorily with RHIC results but not with SPS, although the interesting behavior, at low energy. 

We conclude that HRG model reproduces the experimentally measured fluctuations in (A'+ + K~)/{p + p) and 
{p ~'rP)/{'^^ + ^ ) ratios over the entire range of y/s. In generating this excellent agreement no statistical fitting has 
been performed. Two essential parameters have been merely adjusted with the experimental acceptance in order to 
reproduce the experimental data. 
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SPS 


0.35 


0.1 
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K+ + K-)/{p + p) 


RHIC 


0.35 


0.1 


2.285 




RQMD 


0.8 


0.1 


- 




SPS 


0.55 


0.55 


— 


(p + p)/(7r++7r-) 


RHIC 


0.55 


0.55 


0.636 




HSD 


0.6 


0.55 


- 



TABLE I: Parameters used in HRG model in order to reproduce the experimental and simulation data 

In light of the best reproduetion of dynamical fluetuations that have been measured in various particle ratios over 
a wide range of ^/s, this model, which has been introduced in Ref. [1|, would be a suitable platform to describe other 
collective properties, such as fluctuations in nett charge, and energy, etc. Other implications are also possible. Based 
on the pioneering work about bulk and shear viscosity in hadronic matter [2^, the transport properties turn to be 
accessible by means of of this model. Thus, the comparison with the transport models, like RQMD, UrQMD and 
HSD, etc. shall get more and more creditability. 
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